Cooper pair turbulence in atomic Fermi gases 
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We investigate the stability of spatially uniform solutions for the coUisionless dynamics of a 
fermionic superfluid. We demonstrate that, if the system size is larger than the superfluid coherence 
length, the solution characterized by a periodic in time order parameter is unstable with respect 
to spatial fluctuations. The instability is due to the parametric excitations of pairing modes with 
opposite momenta. The growth of spatial modulations is suppressed by nonlinear effects resulting in 
a state characterized by a random superposition of wave packets of the superfluid order parameter. 
We suggest that this state can be probed by spectroscopic noise measurements. 
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How does a homogeneous interacting many-body sys- 
tem develop spatial modulations as a result of a uniform 
quench? In general, this can happen due to an interplay 
between the extra energy introduced by the quench and 
an intrinsic coupliiig of the degrees of freedom at various 
length scales 0, 0, S 0| ■ nature of the prob- 

lem, the energy distribution of these degrees of freedom is 
initially far from a state of thermodynamic equilibrium. 
This situation is common in a nonlinear medium and can 
be described by the term "wave turbulence" Q. In the 
case when one of the parameters of the medium or an ap- 
plied field is periodically varied in time, wave turbulence 
due to a parametric instability can develop. Examples of 
such a phenomenon include the decay of high-frequency 
electric field into Lagmuir and ion-sound waves in plasma 
01 and the spin-wave instability in an rf-magnetic field 
in dielectric ferromagnets [8,] . A natural question is then 
whether a quench can lead to a parametric excitation of 
spatial modes. 

Oscillating homogeneous fields can be generated by a 
uniform quench without an external drive. One of the 
recent examples where this takes place is a fermionic 
superfluid quenched by a sudden change of the pairing 
strength There are two types of asymptotic states 

the system can reach depending on the strength of the 
initial perturbation lllll2|. The first type has a constant 



value A{t) — As, while the second one is characterized 
by periodic A(i). This analysis does not take into ac- 
count the pair breaking processes and is thus valid only 
for i < Tg, where is the quasiparticle relaxation time. 
Moreover, the emerging asymptotic states are spatially 
uniform: order parameter evolution was obtained as a 
solution of an effectively zero dimensional problem. Re- 
sults of Refs. [m, are thus valid when the system size 
is smaller than the coherence length, L < ^. The descrip- 
tion of the order parameter evolution for i > ^ requires 
an additional investigation. 

In this paper, we perform a stability analysis of the 
solutions for the wave function and order parameter ob- 
tained in Refs. [ll|, [l^l with respect to spatial fluctua- 




tions. We find that the asymptotic states with constant 



FIG. 1; Asymptotic state with the periodic order parameter 
A{t) = As[l -I- gcos(2Ast)] (top left) is unstable with respect 
to the parametric excitation of two modes with opposite mo- 
menta (right panel) . Frequency of each mode is half the oscil- 
lation frequency of the homogeneous order parameter. Initial 
exponential growth of the parametric instability is followed by 
a transient regime after which the condensate reaches a spa- 
tially inhomogeneous post-threshold state (bottom left). The 
value of the order parameter increases from white to black. 



order parameter remain stable, while the state with pe- 
riodic A(i) does not. The physical origin of the instabil- 
ity lies in the possibility of parametric excitation of the 
spatially modulated pairing modes. Fig. [T] In a homo- 
geneous medium, the periodic (in time) order parameter 
can be considered as an energy pump allowing a coupling 
between two spatial modes with opposite momenta, and, 
at the same time, providing enough energy to overcome 
the damping due to the scattering of Cooper pairs. Sub- 
sequent scattering effects of the resonant modes limit the 
initial exponential growth and result in a state with a spa- 
tially inhomogeneous order parameter. We demonstrate 
that as the process of the parametric instability develops, 
the energy of the homogeneous state is transferred into 
that of the pairing wave packets with the typical size of 
the order of the coherence length, ^, signaling the onset 
of the Cooper pair turbulence. Since the amplitudes of 



2 



the wave packets are essentially random, we suggest that 
this state can be experimentally probed by the noise mea- 
surements. 

Consider the state with a periodically varying order 
parameter. Analytically, A{t) is described by the Jacobi 
elliptic function dn Here we assume that the ampli- 
tude of the oscillations is small. This allows us to keep 
only the first two terms in the Fourier series for dn: 



A{t) = As[l + qcos{2Ast)], g<l. 



(1) 



We note that the non-dissipative dynamics within the 
BCS model is described by the Bogoliubov-de Gennes 
equations, which can be cast into the form of equations 
of motion for classical vector variables Sp: Sp = 6p x Sp. 
Here bp = 2(— A(t), 0, Ep) plays the role of an external 
magnetic field. The periodic external field makes it pos- 
sible for the parametric instabilities to develop. Given 
that the asymptotic state dp is robust against homo- 
geneous perturbations 0, [lo( |. it seems natural to look 
for an instability with respect to spatial fluctuations. In 
particular, we will investigate the conditions of the para- 
metric decay of the homogeneous pairing mode ^ with 
an energy 2As into two pairing modes with energies and 
momenta (Li;i,ki) and (a;2,k2), Fig. 1. In a continuous 
medium the energy and momentum have to be conserved: 
ki + k2 = and 2As — uJi + lu2. These type of instabili- 
ties appear in various physical systems, see e.g. Ref. [5j 
for an extensive review. 

To analyze the stability of the asymptotic state ( HI), w e 
employ the Bogoliubov-de Gennes (BdG) equations [13|: 



mp(r, t) = ^Mp(r, t) + A(r, t)vp{r, t), 
Wp(r, t) = -|wp(r, t) + A(r, t)Mp(r, t). 



(2) 



Here ^ ~ —V^ /2m — /i, /i is the chemical potential and 
the order parameter A(r,<:) is 



= g^Wp(r,i)wp(r,t), 



(3) 



where g is the BCS coupling constant. We linearize Eqs. 
^ with respect to the deviations 0p(r,i) and -0p(r,t) 
from the homogeneous solution Up{t) and Vp{t): 
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The order parameter ([3]) also contains a small inhomoge- 
neous part A{r,t) = A{t) + 6A{r,t). Plugging ^ into 
([2]) and Fourier transforming with respect to the spatial 
coordinates, we find 
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The time dependance of Up(t) and Vp(t) (see Ref. fl3|) 
suggests that for linear corrections (|4|) we write 0p (k, t) = 
ap(k,t)e*«p*-F6p(k,i)e-'«p* and Vp(k,i) = apCk,t)e'^''* + 

6p(k,i)e-*«p* with^p = 
/i. We also write 



4 + A2 + 0(g) and £p 



p 

2m 



(5A(k,t) = Ck(i)e 



iAst 



-lAst 



(6) 



where Ck(t) = Cke''(*-*o), Ck{t) = 2ke'''*-*o\ deter- 
mines the growth rate, and to is a time scale when ([T]) 
is reached (in what follows we set to = 0). In the 
expression ^ we have neglected the higher harmonics 
g±i3Ast^ giiSA^t g|-(,^ rpj^jg justified for q <t: I, since 

their inclusion yields higher order in q corrections to the 
growth rate and to the order parameter amplitudes . 
In the linear corrections to the Bogoliubov amplitudes 
(see above) we also keep only the lowest harmonics oj = 
±A„ i.e. ap(k,t) ^ (ai,p(k)e*^=* + a_i,p(k)e-*^^*)e''*, 
6p(k,t) ^ (6i,p(k)e^^^*-H5_i^p(k)e-^^^*)e''* etc. 

Next, we express the Bogoliubov amplitudes in terms 
of Ck and Ck by equating the coefhcients in front of e^'^=* 
in Eq. ([5]). The resulting amplitudes are substituted into 
the self-consistency equation ([3]). One obtains a linear 
system for the variables Ck,cl]^,clj, and Ck. Expressing 
Ck, c*_i^ in terms of Ck and cl^., we derive 

(a;k + i7k)ck + /ikclk = Oi (t^k-«7k)clk + /^kCk = 0, (7) 

where Wk, 7k, and /ik are nonlinear functions of the 
growth rate v. Eqs. ([7]) can be interpreted as the equa- 
tions of motion for a classical field Ck Then, Wk has 
the meaning of the excitation spectrum of this field and 
/ik ^ 0{q) stands for the pumping amplitude, which gives 
rise to a parametric instability. Finally, 7k describes the 
damping of the parametric modes due to the intrinsic 
relaxation processes. 

Nonzero solutions of Eqs. ([7]) for v{k.) exist provided 
Ld'^{v) = h1.[v) — 7k(i^)- Thus, the stability analysis is re- 
duced to the solution of the nonlinear equation for i^(k). 
We have analyzed this equation numerically and present 
the results on Fig. 2. We find that the instability re- 
gion is centered around km ~ l.Qk^ and has a width 
5k « 1.2y/gfc^, where A:^ = As/vf = 1/^ is the coherence 
wave vector. From ^ it follows that for a fixed q the 
parametric growth will be suppressed as soon as the en- 
ergy pumped into the system fully goes into dissipation. 
This condition determines the maximum growth rate v,m 
i.e. 7k (I'm) = h\^{vm)- Our estimate yields Vm ~ 2qAs. 
Lastly, we have also verified that asymptotic states with 
constant order parameter remain stable with respect to 
the spatial fluctuations of the above type. 

The initial growth of the parametric instability (fTTj) 
will be limited by nonlinear effects which lead to the 
transient behavior with subsequent transition into a post- 
threshold state. The latter is defined as a state in which 
Fourier components of the order parameter ^ are time 
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FIG. 2: Region of the parametric instability of the homo- 
geneous A(t) ([T]) with respect to generation of the pairing 
modes with opposite momenta (k, — k). Instabihty growth 
rate is plotted for q — 0.05 in the units of As (see Eq. p)) 
and momentum is in the units of = As/vf- The maximum 
rate is reached at Vm ~ 2gAs. For small q the shape of the 
instability curve is i^{k) — 2Aa(fc — 1.6fc^)^/fc|. 



independent, Ck(t) — Ck and Ck(i) — Ck. Below we focus 
on finding the resulting post-threshold state of the con- 
densate. From the linear analysis we have seen that the 
fastest growing modes are the ones with a certain magni- 
tude of the momentum. Thus, in the BdG equations ([2]) 
among the nonlinear in powers of Ck, Ck terms we keep 
the resonant ones with frequencies lo = ±As and mo- 
menta |k| = |k'| ~ ks, where kg is a new post-threshold 
state momentum to be determined below. The resulting 
set of nonlinear in Ck equations for the order parameter 
amplitudes can be written as Eq. ([7]) with renormalized 
coefficients 

Wk fik = Wk(0) + ^ Tkk'Ick'P, 



k'| = fcs 



Pk^ ft.k(0) + ^ S'kk'Ck'C-k' 

|k'|=fc.. 



(8) 



where Tkk' and 5'kk' are the scattering matrix elements. 
They vary slowly on a scale of fcj and are almost inde- 
pendent of the angle between k and k'. In what follows 
we neglect the fc-dependence in the scattering matrix el- 
ements, Tkk' — T and S'kk' = S. Note that each con- 
tribution in ([8]) is either phase independent or depends 
on a sum of the two phases of Ck and c^j^. This can 
be interpreted as follows. There are two physical pro- 
cesses which limit the parametric excitations: one has to 
do with the reduction of the absolute value of the ampli- 
tudes, while the other is related to the phase decoherence 
of the two pairing modes with opposite momenta. In its 
spirit approximation ([8|) is similar to the mean-field BCS 
model where only diagonal in momentum terms are kept 
in the interaction. The inclusion of off-diagonal terms in 
Eq. ([5]) is expected to cause a broadening of the post- 



threshold state momentum 5k ^ ^/qk^, see Fig. 2. 

To determine the parameters of our post-threshold 
state, we insert Ck — Ickle*""*, and Eqs. ^ for fik and 
Pk into ([7]). The post-threshold state momentum ks is 
determined by the condition that the magnitude of the 
pumping field |Pk| does not exceed the damping 7k for 
any k. As a result we have rik^ = 0. Phase ^'^ — ak+a_k 
and amplitude |c/cj are given by sin^^ = ^kjhk^ and 
|cfe, P = hk, cos 4's/|5'| (the corresponding expressions for 
(5k + ct-k can be derived similarly). We 



|cfcj and *s 
obtain 



ik-r 



A(f, i) + 7^A,c, ^ e- 

\M=ks (9) 

where kg « 1.73A:j, Cg ~ 0.77 and Wg « 0.95 for 
As = 0.1/i. Note that the post-threshold momentum 
kg > k,n ~ 1.6k^, i.e. the energy cascades to smaller 
length scales, as expected of turbulent behavior. 

The individual phases ak and fik cannot be determined 
within the diagonal approximation ([8]). In a continuous 
medium, one can treat them as random variables. For 
the correlators we take (e*"'') 



0, (e'"''! e*"''2 ) = and 
stands for averag- 



^giakig^ka) ^ 4!!-k2e'*% where 
ing over the phase distribution. 

To get further insight into the nature of the post- 
threshold state, consider the following choice for the 
phases ak = — k-ro and ak = ^s/2 — k-ro. It leads 
to a spherically symmetric wave packet, a "bubble" , with 
a periodic amplitude A{t) 

Af^-^\ A , VgAsCsSin(fc,|r-fo|) 

A r,t = Ag + ^ — — — A{t), 

kg\r-ro\ (10) 

A(t) =e^(^+^-^*)+«;,e*(%-^=*). 

In general we obtain a linear combination of these bub- 
bles by writing e*"'' — J f{ro)e~^^''^"d^ro and sim- 
ilarly for e*"''. Then, Eq. ^ can be viewed as a su- 
perposition of wave packets of the form pH)) centered at 
different tq with random amplitudes A{rQ,t). This sug- 
gests that the parametric instability results in a random 
distribution of the wave packets. 

It is instructive to compare (fTO|) with 6A{f,t) at the 
linear stage of the parametric instability, which can be de- 
rived using our result for ^{k) (see Fig. 2) and Eqs. (|6l7p . 
Taking Ck = C'e~*('' '"o+^='^'') and Ck = c^,^, we obtain 



SA{r, t) 



Ce'^'"* cos[As(i - r)] sin(fc„i?)e 



where l{t) 



(«7fc 



(11) 

^ Agt, R = |r — ro|, C is a constant, and 
' LOk)lhk- In deriving Eq. pT|) . we also 



assumed fcmP ^ Agt and replaced a slowly varying func- 
tion Tk — > T. Expression pip describes the initial forma- 
tion of a wave packet pUj) . Note that on a time scale 
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{qAs)~^ at which the order parameter deviation is of the 
order y/qAg, the width of the packet is Ip « H ^/q- 

The observations above help to identify features of 
the post-threshold state ([9|) relevant for the experimen- 
tal verification of our theory. To be more specific, 
let us compute the correlator /C(ri — r2,ii — — 
((5A(fi, ii)(5A(r2, ^2)), where ti and t2 are taken in the 
post-threshold stage. IC{r,t) characterizes the spatial 
and time distribution of the parametric noise in the sys- 
tem. Using the correlators for the random functions e'"'' 
(see above), we obtain from Eq. Q 



JC{ f, t) cx q/S. 



2 sin(/csr) cos A^. {ti 



t2) 



(12) 



where r — jri — r2|- This means that the spatial noise 
spectrum cx J IC{k,uj)du} has a peak at the wave vector 
k — kg and decays as for k ^ kg. 

Formation of an isolated wave packet (fTO| induces 
an oscillating supercurrent js oc V$(r, t), where $(r, t) 
is the phase of the order parameter. Setting tq = 
we find that only the radial component of the cur- 
rent is nonzero. To the lowest order in q, js{r,t) cx 
er^/qcos{IS.st)[ksr cos{ksr) — sin(fcsr)]/(fcsr)^. This im- 
plies a spatial re-distribution of Cooper pairs similar to 
the Friedel oscillations in the density of a degenerate 
Fermi gas induced by a weak scattering potential. 

In our discussion so far we treated the pairing mode 
giving rise to the parametric instability as an inde- 
pendent external field. Inclusion of the feedback on this 
mode as weak turbulence {q <C 1) develops may modify 
the post-threshold state (O . We leave a detailed analysis 
of possible feedback effects for future studies. 

In the post-threshold state ^ the Fourier components 
of the order parameter are Ck,Lj ~ 6{k — ks)6{uj — As). 
Inelastic scattering or thermal effects generally leads to 
a broadeningin the momentum and frequency distribu- 
tions of Ck,uj [al • The latter might cause a damping of the 
temporal oscillations in Eq. ([9]). On a time scale t > 
dissipation due to quasi-particle scattering processes ul- 
timately forces the system to reach an equilibrium state. 
Finally, we comment that in the transient regime leading 
to an asymptotic state with constant A{t) = As the order 
parameter is (A(t) - A^) cx cos{2Ast)/ ^/i (cf. (P)) [l3|. 
Oscillatory behavior suggests that this asymptotic state 
might never be attained owing to the development of the 
parametric instability of the type considered above. 

In conclusion, we have investigated the stability of the 
nonequilibrium asymptotic states of a fermionic super- 
fluid, which can be generated e.g. by a uniform quench 



of the pairing strength. We have demonstrated that in 
a system of size L larger than the coherence length ^ 
the asymptotic state (P) with periodic in time order pa- 
rameter is unstable with respect to spatial fluctuations. 
The instability is due to the parametric excitation of two 
pairing modes with opposite momenta. The initial expo- 
nential growth of deviations from the homogeneous state 
(fTTj) is suppressed by nonlinear effects eventually leading 
to a spatially nonuniform post-threshold state described 
by Eq. This state can be interpreted as a superpo- 
sition of bubbles of the superfluid order parameter (jlOp 
with random amplitudes. The parametric instability of 
the uniform oscillations can be experimentally probed by 
spectroscopic noise measurements. 
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